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Abstract 

We have derived a model of the Kuiper belt luminosity function exhibited by a 
broken power-law size distribution. This model allows direct comparison of the 
observed luminosity function to the underlying size distribution. We discuss the 
importance of the radial distribution model in determining the break diameter. We 
determine a best-fit break-diameter of the Kuiper belt size-distribution of 30 < 
D}j < 90 km via a maximum-likelihood fit of our model to the observed luminosity 
function. We also confirm that the observed luminosity function for m{R) ~ 21 — 28 
is consistent with a broken power-law size distribution, and exhibits a break at 
m{R) = 26.01?:^. 
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1 Introduction 



The size distribution (SD) of a belt of planetesimals can provide information 
on the physical and dynamical properties of those objects. Colhsions between 
belt members causes the SD to evolve over time. Low velocity collisions can 
cause objects to accrete into larger bodies, while high velocity collisions smash 
bodies into smaller pieces. A generic outcome of coUisional evolution simula- 
tions of a belt of planetesimals is a power-law SD with a steep slope for large 
objects, breaking to a shallower slope at some object diameter. Waves or other 
shapes can be induce d in a size distribution, but the general behavior is that of 



a broken power-law (jO'Brien and Greenberg 



20031 ). The location of the tran- 



sition between the steep and shallow SD slopes (break), as well as the small 
and large object slopes provides information on the strengths of the collid- 
ing bod ies, and constrains time-scale s of accreti o n and total mass in the belt 



region fISafronov anc 


2001; 


Kenvon. 


20oJ: 



1969 



1 r 



Dohnanvi 



Kenyon and Bromleyl . 



2004 



119691: iKenvon 



Bottke et al. 



and 



Bromley 



20051 ) 



The SD of the Kuiper belt has been inferred from a measurem ent of the 
belt's luminosity function (LF) (see for example 



Fraser et al. 



2008 



(F08) and 



references therein). FOB show that for the largest objects in the belt (diameter 
D > 100km), the SD is well represented by a power-law N{D) oc D^'^ with a 
slope g ~ 4.3. 



Kuiper belt surveys designed to measure the Kuiper belt object (KBO) faint- 
end LF (m(i?) > 26), have found eviden c e to suggest that the LF deviates 



from a power-law (IBernstein et al. 



20M). 



Bernstein et al. 



(12004h found that 



a distribution with a steep slope for bright objects, that "rolled over" to a 
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shallow slope for faint objects was necessary to describe the substantial dirth 
of objects with m{R) > 27. They concluded that the underlying KBO SD 
transitions ("breaks") to shallower slopes as predicted f rom coUisional evolu- 



tion m odels. Using archival Subaru Suprimecam images, 



Fuentes and Holman 



( I2OO8I ) provided a substantial increase in the sample of faint {m{R) ~ 26) 
KBOs and concluded that the luminosity function breaks from a steep to 
shallow slope at m{R) ~ 24.3. 



Bernstein et al 



tQ04 ] and 



Fuentes and HolmanI (l2008l ) use analytic functions 



to represent the shape of the observed LF. These functional forms were not 
derived from physical or theoretical properties, making infering the underlying 
KBO SD based on the reported LF difficult. We present here a new functional 
form of the LF that results from the observation of a belt of planetesimals 
exhibiting a broken SD. This function provides a direct interpretation of the 
observed LF in terms of the size and radial distributions. 



In Section 2, we derive the new functional form. We fit the function to the 
available observations of the Kuiper belt in Section 3, describe what future 
efforts are necessary to constrain the SD in Section 4, briefly discuss impli- 
cations of the Kuiper belt SD in Section 5, and make concluding remarks in 
Section 6.. 



2 Broken Power-Law Derivation 



To determine the shape of the LF exhibited from a belt of planetesimals with 
a given SD, we flrst consider a number of objects in a sample represented by 
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N c:^A R{r-)S{D)dDdr 



(1) 



where r is heliocentric distance, D is object diameter, and R{r) and 5* {D) 
are the radial and SDs of those objects (number of objects with distance r to 
r + dr oi diameter D to D + dD). A is a normalization constant to present 
the equation in number of objects per square degree. Here we have implicitly 
assumed that the radial and size distributions are separable. 

The magnitude of brightness of an object at a heliocentric distance r in AU, 
with diameter D in km is well approximated as 



m 



ir + 2.51ogio(r^I^-') (2) 



where, for KBOs we have applied the simplification, A ~ r. Setting K = 
8.4, corresponds to KBOs with 6% albedos when observed at opposition 



(Gladman et al. 



200ll ): implicitly assuming that all objects have similar albe- 
dos. The largest object observable at a given distance, r in a survey that is 
sensitive to source magnitudes between m^m and m^ax has diameter D^ax = 
7-2]^g(m,mi„-ft:)/5_ gjxxiilarly, the smallest observable object has diameter Dmin = 

^2'^Q{mmax—K)/5 

The cumulative LF of a belt of planetesimals as observed by a particular 
survey, can be derived from Eq. 1 by using the integration limits Dmin and 
Dfnax- If we assume S{D) oc ie. we assume a SD is a power-law with 

'slope' q, Eq. 1 becomes 



5 



N (< nir, 



A 



dr 



10" 



-{l-q) 



il-q) 



(3) 



where 7^ 1. Setting a = ai = A / ^-^^dr^ and restricting to the case, 
10 1-9 ^ 10 , then the cumulative LF can be represented by 



N {rrimax) = Oi {rrimax) 10° 



(4) 



We consider the case where the SD breaks to a different slope at a particular 
break size Di,. That is, 



oc 



for D > A 
for D < Db. 



(5) 



If ri and r2 define the inner and outer edges of the belt, then a survey where 
Dmin = rilO^"*"*"^"-^)/^ > Db, the observed LF is given by Eq. 4 where the 
slope a — ai — If however, the survey is sensitive to objects smaller than 
Db out to some distance Vb, then Eq. 1 becomes 



N 



A R{r) (5 D-i'dD+ D-^'dD 

Jri JDmir,. JDh 



^max 



dr+A / R{r) / 



D-'^'dDdr 



(6) 



where /3 is a scaling factor to ensure that S{D) is piecewise continuous at Db- 

Using Eq. 2, the maximum distance a survey is sensitive to the break size 
Db, is given by n = L)|lO("*'"»---^)/^°. 



The constant /3 can be solved for by requiring that the number of objects 
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slightly larger than Db be equal to the number of objects slightly smaller than 
Db. That is 



JDh-S JDb 



(7) 



which gives 



(8) 



Restricting to the case where 10^^(i-«i) < 10^"T^(i-5i) and lo^^^(i-?i) < 
10 's™'""" (^~'J2) ^ then the cumulative LF is given by 



where 



Qo — 1 qi — 1 

N (< m^ax) - 6ilO— + 6210—"^— + 63 {rrimax) (9) 



and 



Jri q2 — ± 



,2(l-qi) 



62 = A / i?(r) drlO— 

^r;, Ol - 1 



.£)l/2lo('"r"ax-if)/10 
6 



g2 - qi 



(gi - 1) (g2 - 1) 



&3 = ^ / " R{r)drD^- 

where qi 1 and 52 7^ 1, and A is chosen such that N is the number of objects 
per square-degree of the sky. 



If a survey is deep enough to see objects with diameters D < Db for all r, 
then the cumulative LF then becomes 
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where 



and 



N{< rriraax) ^ CilQ— + C2 (10) 



rr2 y-2(l-q2) -K(qy-l) 

ci^A R{r)- -Dl'-'^'drlO^^ 



Jri (1 - q2) (1 - qi) 

The cumulative LF of a belt of planetesimals with a two sloped power-law SD 
is a power-law with slope cti for magnitudes brightward of mi. The roll-over in 
the LF begins at the brightest magnitude at which objects smaller than can 
be detected at the inner edge of the belt, ri. Thus, mi = K + 2.5log (rfD^'^^. 
The LF becomes a power-law again at the magnitude m2 faint-ward of which, 
objects larger than the break-size are no longer observed at the outer edge of 
the belt, r2. Thus, m2 = -f^ + 2.5log ^r2-D^^). For magnitudes mi < m < m2, 
the LF is the sum of nearby objects with D < Db, and distant objects with 
D > Df,. Hence both the radial and size distributions become important as 
the shape of the LF is determined by the coupling of the two. The functional 
form of this LF is 



N(<m) = A 



ailO"!"* if m < mi 

6il0"2™ + 6210^1™ + 63 (m) ifmi<m<m2 (^1) 

Ci + C2lO"2'" ifm>m2. 



Comparing to Eq. 9, we see that the SD slopes qi and q2 are related to the 

LF slopes cti and 0:2 by qi = 5ai + 1 and q2 = 5a2 + 1. 

The power-law coefficients, a^, 6^, and q, of this function depend on the radial 
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distribution. For the case where a break is observed, knowledge of the true 
radial distribution is important as the shape of the LF is determined from 
the coupling of the radial and size distributions. For the case where the break 
is not seen, the radial distribution is not important as it only affects the LF 
normalization. The size distribution slopes qi and q2 can be inferred from the 
LF slopes «! and 02 if the LF has been observed sufficiently brightward and 
faint-ward of the roll-over. 



3 Model fit to the Observations 



We present here a fit of our model (Eq. 11) to the observed KBO LF. 



Eq. 1 can be evaluated analytically if R{r) oc r where c is the radial dis- 
tribution 'slope'. While t his is likely not the true KB O radial distribution. 



Trujillo and Brown 



200 If ) and 



Kavelaars et al 



(I2OO8I ) have shown that the 



radial distribution is s harply decreasing beyond r ~ 40 AU. A visual examina- 
tion of Figure. 2 from lxrujillo and Brown ( 2001 ) reveals that R{r) = (^-^^ 



is a good representation of the sharp density fall-off. Such a representation 
is sufficient to show the importance of accurately modeling the radial distri- 
bution when determining the size distribution. The full functional form using 
this radial distribution is presented in the appendix. 



Eraser et al. 



(I2OO8I ) present a fit of the LF for data brightward of m{R) = 26. 
They find that the LF of the Kuiper belt is well described by a single power-law 
with slope a = 0.65 ±0.05. To ensure reliability of the fit, they only considered 
surveys with measured detection efficiencies as a function of magnitude for 
each field in the survey. In fitting Eq. 11 to the observed LF, we adopted the 
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same practices as iFraser et al\ (120081 ) listed below: 



tion, we included the survey presented by 



/ Fraser et al. 


( 


2008 


).I 


n addi 


Bernstein et al 


• ( 


2004. 


2006) 



Note: we did not included the survey presented by 



Fuentes and Holman 



( I2OO8I ) as the majority of analysis for this manuscript was complete be- 
fore the author was aware of that work. This manuscript focuses on the 
functional form of the LF, and these data, which may improve the con- 
straint on the function parameters, do not change our conclusion that a 
full understanding of the size distribution requires a rigorous treatment 
of the radial distribution as well. 
(2) We offset all magnitudes to R-band using the average KBO colours pre- 



sented in 



Fraser et al. 



mm . 

(3) We cull from each survey all sources faintward of the 50% detection ef- 
ficiency of that survey, and set the detection efficiency to zero faintward 



'^^drn^^ using a maximum likelihood 



of that point. 

(4) We fit the differential LF, E(m) 
technique. 

(5) We adopted a form of the likelihood equation as derived by 



LoreddfeoO^ 



and extended to account 



tions, Ch fseelFraser et al. 



' or sky density variations, A^, and color varia- 



20081 . for details). We adopted the same prior 
ranges of the colour offsets an d range of normalization parameters as pre- 
sented in 



Fraser et al. 



(120081 ) (±0.2 for the colour parameters, and ±0.4 



for the normalization parameters). 

To evaluate Eq. 11, we require a set of radial distribution parameters, ri, r2, 7 = 
(1 — c)/10. We set the Kuiper belt inner edge as ri = 35 AU. For a baseline 
model, we consider an outer edge r2 = 60 AU with the radial fall-off parameter. 
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c = 10 (see previous discussion). This provides a reasonable representation of 
the steep fall-off in the Kuiper belt radial distribution and typical outer limit 
to the objects detected in the F08b sample. We also consider two test cases; 
in model 2 we set r2 = 100 AU, and in model 3, we set the radi al fall-off 



parameter, c, to that of the minimum-mass solar nebula, c = 3/2 (IHayashi 



19811 ) corresponding to 7 = —0.05. Models 2 and 3 were chosen to evaluate the 
effects of the radial distribution parameters on the fit results. We speculate 
that the value of r2 for any reasonable choice, will have no strong effect on the 
fit parameters as most objects are detected at distances less than ~ 45 AU. 

Using our maximum likelihood technique, we determined the best-fit LF pa- 
rameters of Eq. 11 using the baseline model, as well as models 2 and 3. These 
are presented in Table 1. The observations were best-fit with ai = 0.69 ±0.07, 
Db = 371^0 km, a2 = — 0.4l3^ for our baseline model. 

We found that the best-fit ai and Db, Db = 371*110 km and ai = 0.69±0.07, are 
not dependent on the choice of radial distribution models. As demonstrated 
in Eqs. 3 and 9, the radial distribution has no effect on the apparent bright 
object slope and the break diameter depends only on the choice ri. Thus, 
because the same data was always used, and because we set ri = 35 AU for 
all choices of radial distribution, the best-fit ai and Db should be the same for 
all radial distributions considered. 

We found very different behaviors for the faint-end slope from our two choices 
of 7. The best-fit 0:2 was found to be — O.4I3'®, — 0.4l3'^, and —00 for the 
baseline model (n, r2, 7) = (35, 60, -0.9), Model 2 (n, r2, 7) = (35, 100, -0.9), 
and Model 3 (ri,r2,7) = (35, 100, —0.05) respectively. 

With 7 = —0.05 (c = 3/2), we found that the observations were best-fit with 
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a SD in which there is a complete absence of objects, {a2 — — oo) smaller than 
break diameter, Db ~ 40 km. When 7 = —0.9 (c = 10), the best-fit faint end 
slope was a — —0.4. These two behaviors are expected; a sharp roll-over in 
the LF can be the result of a very sharp roll-over in the SD with a gradual 
fall-off in the radial distribution, or a less sharp roll-over in the SD and a steep 
fall-off in the radial distribution. Setting c — 3/2 does not correctly model the 
very sharp radial fall-off in the Kuiper belt and thus it is unlikely that there is 
a complete absence of objects smaller than ~ D^. We discard model 3 as this 
is an unlikely model of the radial and size distributions. We do note however, 
that no KBOs smaller than D ~ 30 km have yet been detected. 

We found no difference in the best-fit parameters when r2 was varied (see the 
basehne model and model 2). This demonstrates the expected lack of influence 
of this parameter on the fits. 

A series of Monte-Carlo simulations were performed in which a sample of 
detections representative of the surveys and detections in the F08 sample were 
generated from a model Kuiper belt. Each simulated data-set was fit using the 
maximum likelihood technique, and on average, the fitting routine reproduced 
the input model parameters. The same behavior was found for the fits of the 
simulated observations, as for the real observations. That is, the best-fit faint- 
end slope, a2, depended on our choice of 7, and the best-fit break diameter, 
Dfe, depended on our choice of r^. The reported best-fit values were correct 
however, if the chosen 7 and ri were within a factor of 2 the correct value. 
Thus we feel that our likelihood method presented here produces satisfactory 
results to within the accuracy of the fit to the data currently available. 

All minimizations have the same maximum likelihood value to within a few 
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percent. Thus there is no preference, in a statistical sense, to choose one 
model over another; each is equally sufficient to describe the observations. 
Thus, we take take the baseline model, with (ri,r2,7) = (35,60,-0.9) and 
(logy4, ai, 0:2, -Dfe) = (22.82,0.69,-0.4,37) as the best model, as the baseline 
model radial distribution is the most similar of the three considered, to that 
observed for the Kuiper belt. 

The differential LFs (derivative of Eq. 11 wrt. m) using the best-fit parameters 
for all choices in radial distribution are presented in Fig. 1. To produce the 
histogram, object magnitudes and the areas of each survey have been adjusted 
using the best-fit Ck and Ak parameters (see F08 for details). These parameter 
values were determined by maximizing the likelihood while holding the LF 
parameters to their best-fit values, and are presented in Table 2. The histogram 
is presented in order to visualize the observations, but was never used in 
the fitting procedure. As can be seen, a very sharp roll-over is apparent for 
magnitudes faint-ward of m{R) ~ 26. This plot provides visual con&mation 
that the fit describes the data reasonably well. 

The likelihood contours of the best-ffi LF for our choice in baseline radial 
distribution is presented in Fig. 2. As can be seen, ai and A and 02 and Df, 
are highly correlated. The uncertainties of the best-fit parameters in Table 1 
are the extrema of the 67% bayesian credible regions. 



The best-fit large object sl ope and break diam eter, qi = 4.45 ±0.35 and Db 



371^0, are consistent with (IFraser et al. 



20081 ) who found that q = 4.25 ± 0.75 



for objects larger than D ~ 50 km. 
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4 Sources of Uncertainty 



The best-fit LF we pr esent assumes that the lack of fai nt KBOs observed by 



Bernstein et al 



fl2004h and 



Fuentes and Holman 



(I2OO8I ) is caused by a break 



in the SD, and suggests that the break observed must occur at diameter 27 < 
Dh < 87. This resuh however assumes that an R-albedo of 6% is typical of 
KBOs. Recent observations suggest that this is true for object s D < 500 km , 



but 


ob 


2005 


1 ^ 



Stansberry et al. 



20071 ). 



I I I 

Fraser et all (i2008i) has shown however, that 



the trend of lower albedos for smaller objects would produce a steeper LF 
than otherwise would be observed for constant albedos. Thus, albedo effects 



are likely not the source of the observed break. In addition. iFraser et ali (12008 ) 



found that a variation of albedo with distance has no effect on the inferred 
LF slopes, and only the LF normalization. 



We also assume that the size distribution is the same at all distances. If this 
is not true, it would have similar affects as assuming constant albedos for all 
sizes. That is, the slope inferred under the assumption that the size distribu- 
tion is distance dependent could be steeper or shallower than the true slope 
depending on the trend of size with distance. This effect cannot be tested 
without substantially more observations than currently available. This effect 
however, is likely negligible, as KBOs have excited orbits. Thus interactions 
(collisions) between objects occur over a range of many tens of AU which would 
cause the size distribution to be similar over the full range of the Kuiper belt. 

The LF also assumes that the break in the SD is sharp. That is, there is 
no size range over which the SD slope transitions from qi to q2- While this 
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is certainly not true, simulations of accretion and coUisional processing in 
a belt of planetesimals show that the transition occurs over a small range 
of sizes c ompared to th e range over which the SD exhibits a power-law like 



behavior (IKenyon 



20021 ). Additionally, the lack of surveys sensitive to KBOs 
with m{R) > 27 from which the roll-over shape is constrained implies a more 
complicated description is not warranted. More survey data is needed that 
covers a few square degrees in the range m{R) ~ 26—28 to accurately constrain 
the break magnitude; a large number (A^ ~ 20) of detected objects with a 
range of magnitudes beyond the break are necessary to constrain the faint end 
slope, without which the break diameter determination remains uncertain. 



The SD inferred from the LF is still model dependent, and determining the 
structure of the SD reliably from the LF requires improved knowledge of the 
properties of the Kuiper belt (albedo, radial distribution, etc.). A large source 
of uncertainty stems from the current lack of knowledge of the KBO radial 
distribution. Each survey considered in this analysis sampled the Kuiper belt 
at different locations on the sky. Variations in the sky density between sur- 
veys, requires tr eating sky density parameters as nuisance parameters for each 



individual field (IFraser et al. 



20081 ) . This causes the range of parameter values 
consistent with the observations to be grossly enlarged. An accurate Kuiper 
belt model which provides the instantaneous radial distribution as a function 
of latitude and longitude is needed such that density estimates and the coeffi- 
cients of Eq. 11 can be calculated as a function of latitude and longitude. This 
would remove the need for treating sky densities as nuisance parameters, and 
would greatly reduce the uncertainty in the inferred SD. 
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5 Discussion 



Our results have demonstrated that the size distribution of the Kuiper belt can 
be described by three size distribution parameters, gi, q2, and Df,, and that, to 
accurately measure the break-slope q2 requires that the radial distribution be 
accounted for in a way which removes the ambiguity between source brightness 
and source size. We have found that gi ~ 4, g2 ~ — 1, and Df, ~ 40 km. 



As in pr evious works, we have found that th e obse r ved large object slope 



1 

Q^ ~ 4 q 


— - "w^^^^, 
ruiillo et al. 




2001 


Petit et al. 


2006; 


Eraser et al. 



Gladman et al 



2001 



Bernstein et al. 



2004 



20081 ) , is incon sistent with the q ~ 3.5 slope ex- 



pected from a steady-state collisional cascade ( IDohnanyi 



2003 



Pan and Sari 



1969 



O'Brien and Greenbergl . 



20051 ). Rather, the slope is consistent with the large ob 



iect slope produced from accretionary 



Kenyon and Bromleyi . 



2001 



Kenyon 



proce sses (ISafronov and Zvjaginal . 



1969 



2OO2I) 



Various accretion and collisional models exist that can reproduce some of the 
observed features of the size distribution, ie. the correct small object slope, 
or break diameter. But t o this date, none ca. n simultaneously reproduce gi. 



g2, and D^. For instance. 



Pan and Sari 



(120051 ) present a collisional model in 



which the accretion size distribution breaks to one in a state of collisional equi- 
librium. They demonstrate that a break diameter consistent with our results 
could occur if collisional disruption occurred for timescales ~ 1 — 4 Gyr. The 
equilibrium slope they considered (g2 ~ 3) however, is i nconsistent w i th the 



Pan and Saril (120051 ) 



small object slope of our baseline model. The assumption 
make, that the size distribution immediately breaks from an accretionary slope 
to collisional equilibrium is likely incorrect, as there is likely a range of sizes 
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over which disruptions by the more numerous smaller objects have significantly 
depleted the population, but which have not been replenished with fragments 
from larger object disru ptions. This transition rang e is seen in more co mpli- 



cated accretion models (IKenyon and Bromleyl . 



2001 



Bottke et al. 



20051). 



In early times of the accretion simulations of 



Kenyon and Bromleyl (120011 ) 



KenyonI (120021 ) . steep large object slopes, gi ~ 4 are produced. For times longer 
than ~ 70 Myr (~ 1 Gyr for weak bodies), the slope rapidly evolves to a shal- 
lower slope inconsistent with our findings. This suggests that accretion in the 
primordial Kuiper belt was a short-lived process which must have been 'turned 
off' by some process - likely the dynamical excitation of the belt by some 
scattering pr ocess, such as gravitational stirring from a migrating or scatter- 



ing N eptune (IMalhotra 



1993 



20031 ) or a rogue planet (iGladman and Chanl . 



Thommes et al. 



2002 



Levison and Morbidelli 



20061 ) . These simulations also 



produce a break diameter, D}^ ^ 10 
results. 



sm - only moderately consistent with our 
Kenyon and Bromleyl (l2004j ) found that, with gravitational stirring 
from Neptune, after 4.5 Gyr, 'weak' KBOs - that is, bodies with collisional 
strengths lower than expected from icy planetesimals - exhibited a roll-over at 
break diameters as large as Di, ~ 30 km, consistent with that found here. In 
addition, there simulations produced small object slopes consistent with our 
findings, but produced a large object slope which is too shallow compared to 
the observed slope. 



he m odels studied by 



Kenyon and Bromleyl (120011 ): 



KenvonI (l2002h : 



Kenyon and Bromley 



( l2004j ) have suggest that KBOs are effectively strengthless, gravitationally 
bound rubble-piles. These models however, cannot simultaneously produce 
both the large object and small object slopes; either the evolution is short 



enough such that the steep accretion slope is maintained, or the evolution is 
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long enough such that a large enough break, and shallow enough small object 
slope are produced. This is not a surprise however, as it is apparent that the 
Kuiper belt history is not as quiescent as that considered in these models. The 
highly excited KBO orbital distribution cannot be produced by in-situ forma- 
tion. Some large-scale scattering event is needed. Thus, the correct dynamical 
history must be determined before reliable conclusions can be drawn about 
the accretion/disruption history, and strengths of KBOs. 



6 Conclusions 



We have developed a model of an LF exhibited from a planetesimal belt with a 
broken power-law SD. This model allows a direct interpretation of the observed 
LF in terms of the underlying SD assuming a broken power-law shape. We 
fit this model to the o bserved Kuiper belt LF and the break first observed by 



Bernstein et al 



(I2OO4J ). The best-fit parameters are presented in Table 1. We 
find that the break in the SD must occur at diameters 27 — 87 km {m{R) ~ 
25.8) for reasonable choices of the radial distribution. We have demonstrated 
that the inference of the SD from the LF requires an accurate model of the 
radial distribution. Additional detections of KBOs faintward of m{R) = 26, 
and a proper treatment of the KBO orbital distribution are required to provide 
an accurate measure of the break size and faint-end slope. 
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8 Appendix 

As an approximation to the true Kuiper belt radial distribution, we assume 
that R{r) oc and ri < r < r2. Making the substitution 7 = (1 — c)/10, 
and collecting like terms, the coefficients of Eq. 11 become 



10(7-«l) _ ^10(7-0!) 

a, = — 10-"^, (12) 

50a 7 -a ' ^ ^ 



D 



5(a2-ai) 



^ -A^^^ -rf (™h0-^'^^ (13) 

50a2 (7 - a2) 



10(7-01) 

h ~ A-^ ^10--^^, (14) 

50q;i(7 — ai) 



63(m) ^ A ^ K^^-^^) Rb ^Q7(™-i^) (15) 

^ (7 - ai) (7 - 0:2) 5O7 ^ ^ 



107 107 , 



and 
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10(7-02) _ 50(7-02) 

'2 'i n''''("2-«i)i n-o2X 



where 7 7^ ai, 7 q;2, cti ^ 0, and a2 ^ 0, or equivalently 3 — 2gi — c 0, 
3 — 2^2 — c ^ 0, 1 — gi ^ 0, and 1 — ^2 7^ and we have dropped the neghgible 
terms from the coefficients. 
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Table 1 

Maximum Likelihood Parameters. The bright and faint object slopes ai, a2, the 
break diameter Di,, and the normalization constant A were fit by the maximum 
likelihood routine. Belt edges ri and r2 as well as the radial fall-off parameter 7 are 
fixed parameters not fit by the maximum likelihood routine. Fit 1 represents our 
baseline model, as the chosen radial distribution best matches that of the Kuiper 
belt. Fits 2 and 3 have different values of r2 and 7 respectively, chosen to determine 
the sensitivity of different choices in radial distribution models on the measured LF 
parameters. Note: 7 = —0.05 produces physical implausible results, ie. a2 = —00. 

Parameter Baseline Model Model 2 Model 3 

n 35.0 35.0 35.0 

r2 60.0 100.0 100.0 

7 -0.9 -0.9 -0.05 

log A 22.8 ±0.6 22.8 ±0.6 9.3 ±0.6 

ai 0.69 ±0.07 0.69 ± 0.07 0.69 ±0.07 

a2 -0.4l^-« -OAtl-^ -00 

D,{km) 371™ 37t™ 37^™ 
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Field \ogAk Ck 



Source 



UN 



22.6 0.0 



MEGA 23.0 -0.15 



Eraser et al. f2008l 
Eraser et al. (2008) 



CTIOOl 23.2 0.1 



Eraser et al. (2008 



CTIO02 23.1 0.05 



Eraser et al. f2008 



SSU 



23.0 0.0 



Petit et al. f2006 



SSN 



GOl 



23.1 0.03 



22.9 -0.01 



Petit et al. (2006 



Gladman et al. (2001) 



AE 



23.0 0.0 



Allen et al. (2002 



AKL 



Bern 



23.2 0.0 



Allen et al. (2002 



TEIG 22.8 0.0 Truiillo et al. (2001) (see Eraser et al. (2008) ) 



TE2G 22.8 0.0 Truiillo et al. (2001) (see Eraser et al. (2008 



TE3G 22.7 0.0 Truiillo et al. (2001) (see Eraser et al. (2008) ) 



TE5M 22.8 0.2 Truiillo et al. (2001) (see Eraser et al. (2008) ) 



22.7 -0.2 



KeckV 22.6 0.2 



JLT98 22.5 -0.2 



Bernstein et al. (2004) 



Chiang and Brown (1999) 



Jewitt et al. (1998) 
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G98UH 22.4 -0.1 



Gladman et al. fl998 



Gladman et al. fl998 



G98Palomar 23.2 -0.03 



TEIM 22.8 0.13 Truiillo et al. f2001) fsee Eraser et al. f2008') ) 



TE3M 22.9 0.19 Truiillo et al. f20011 fsee Eraser et al. f2008) 1 



TE4G 22.8 0.0 Truiillo et al. f20011 fsee Eraser et al. f2008 



TE4M 22.7 0.05 Truiillo et al. f20011 fsee Eraser et al. f2008) 1 



AB 



AD 



AE 



AG 



AJ 



22.9 -0.11 



22.8 0.0 



22.9 0.0 



22.5 0.0 



22.8 0.08 



Allen et al. f2002 



Allen et al. (^002) 
Allen et al. (^002) 
Allen et al. (2002) 



Allen et al. f2002 



Table 2: Best-fit normalization parameters and 
colour off-se t para meters Ck for each field as defined in 



Eraser et al 



(120081 ) . These are treated as nuisance param- 
eters and are marginalized when determining the best- 
fit LE (see Section 3). Eor this table, and Eig. 1, these 
best-fit parameter values are determined from the maxi- 
mum likelihood routine when the LE parameters are set 
to their best-fit values (Eit 1 from Table. 1). 
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Fig. Captions 



Fig. 1 Differenti a 



Bernstein et al. 



surf ace density histogram of data from the F08 sample and 



(120041 ) using 0.8 mag bin-widths. Objects magnitudes have 



been shifted using the best fit colour offsets Ck, and observed sky densities 
have been adjusted using the normalization values taken from Table 2. 
Das hed straight line: be st fit power-law with a = 0.65 and m/j = 23.42 taken 
from lFraser et ali (120081) . Solid curve: Model 1 (basehne model), dotted curve: 
Model 2, and dashed curve: Model 3. The baseline model has the radial distri- 
bution which best represents that of the Kuiper belt. Models 2 and 3 test the 
sensitivity of the fit to different radial distributions. Parameters listed in Ta- 
ble 1. As can be seen, a strong break is seen in the LF, as a deviation from the 
power-law behavior faint-ward of m{R) ~ 25.8. The shape of the binned LF is 
highly sensitive to the choice of bin boundaries and readers are cautioned from 
drawing conclusions based on such representations. We present the observa- 
tions in this form only to demonstrate visually the generally correct shape of 
the LF as derived from our model SD. 

Fig. 2 1,2, and 3-sigma likelihood contours for pairs of bright and faint ob- 
ject slopes ai, and 0^2 (or equivalently large and small object size distribution 
slopes qi and ^2), break diameter (assuming 6% albedos), and normaliza- 
tion parameter log A. Contours for parameter pairs were generated with the 
other pair of parameters set to their best-fit values. 
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Fig. 1. 
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